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Ein  Fliissigkeitsnutationsdampfer  fur  den  Satellit 

EQUATOR-S 


Bernd  Hausler*,  Neubiberg,  Germany 
Werner  Eideb,  Neubiberg,  Germany 

February  2,  2001 


Zusammenfassung 

Zur  Auslegung  eines  Nutationsdampfer  fur  den  EQUATOR-S-Satelliten  wurde  ein  Nu- 
tationsdampfersystem  mit  zwei  gleichen  symmetrisch  angeordneten  Fliissigkeitsdampfern 
betrachtet.  Ein  solcher  Dampfer  besteht  aus  zwei  Rohren  mit  kreisformigem  Querschnitt 
die  durch  zwei  kurze  Rohre  mit  wesentlich  grofierem  kreisformigem  Querschnitt  verbunden 
sind.  Zu  den  vorgegebenen  Tragheitsmomenten  des  Satelliten  sowie  dessen  Spinfrequenz 
wurden  dazu  unter  Beriicksichtigung  gewisser  Einschrankungen  beziiglich  der  Geometrie 
der  Dampfer  als  auch  der  Wahl  der  Fliissigkeit  die  Dampferparameter  so  berechnet,  dab 
eine  gute  Nutationsdampfung  mit  eingeklappten  und  ausgeklappten  Booms  gewahrleistet 
war.  Die  Dampfer  wurden  nach  den  so  berechneten  Parametern  gebaut  und  an  einem 
luftgelagerten  Pendel  getestet.  Hierzu  wurde  auch  die  theoretische  Dampfung  am  Pendel 
gerechnet  und  mit  dem  Experiment  verglichen.  Die  Ergebnisse  von  Theorie  und  Ex¬ 
periment  am  Pendel  stimmte  dabei  gut  uberein.  Eine  Aufzeichnung  des  abklingenden 
Nutationswinkels  an  Bord  von  EQUATOR-S  zeigte  ebenfalls  gute  Ubereinstimmung  mit 
der  Theorie.  Um  diese  gute  Ubereinstimmung  zu  erzielen,  war  es  jedoch,  im  Gegen- 
satz  zu  fruheren  theoretischen  Untersuchungen,  notwendig  eine  Korrektur  der  bewegten 
Fliissigkeitsmasse  im  Dampfer  vorzunehmen. 


*Univ.Prof.Dr.rer.nat,  Universitat  der  Bundeswehr  Miinchen,  Institut  fur  Raumfahrttechnik,  Werner- 
Heisenberg-Weg  39,  85577  Neubiberg 

tDr.-Ing.,  Universitat  der  Bundeswehr  Miinchen,  Institut  fur  Raumfahrttechnik,  Werner-Heisenberg- 
Weg  39,  85577  Neubiberg 
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1  Einleitung 


Am  2.  Dezember  1997  wurde  der  vom  Max-Planck-Institut  in  Garching  gebaute  spin- 
stabilisierte  Satellit  EQUATOR-S  in  eine  Erdumlaufbahn  geschossen  um  Plasma-  und 
Feldmessungen  in  der  Magnetosphere  im  Rahmen  des  ”  International  Solar  -  Terrestial 
Physics  Program”  (ISTP )  durchzufiihren.  Die  Orbitdaten  nach  Ziinden  des  Kickmotors  in 
einem  Ariane  4  Transferorbit  waren:  497  km  Perigaumshohe,  67  275  km  Apogaumshohe 
(11,5  Re  geozentrischer  Abstand)  3,9°  Inklination,  22,3  Stunden  Umlaufdauer.  Damit 
wurde  die  Magnetosphere  der  Erde  in  niedriger  magnetischer  Breite,  deren  Grenzschicht, 
die  dquatoriale  Ringstromungsregion  und  die  erdnahe  aquatoriale  Plasmaschicht  erforscht 
[1.  2]. 

Wie  bei  alien  spinstabilisierten  Satelliten  mufite  auch  bei  EQUATOR-S  gewehrleistet 
sein,  daB  durch  auBere  Storungen  auftretende  Nutationsbewegungen  des  Satelliten  aus- 
gedempft  werden.  Dies  kann  auf  verschiedene  Weise  geschehen:  aktiv  durch  Kaltgasdiisen 
oder  passiv  durch  Dempfungseinrichtungen.  Solche  passiven  Dempfungseinrichtungen 
sind  im  Satelliten  eingebaute  gedempfte  Schwingungssysteme,  die  beziiglich  ihrer  Eigen- 
frequenz  und  Dampfung  so  ausgelegt  werden,  dafi  im  gekoppelten  System  Satellit- 
Dhmper  die  Nutationsschwingung  moglichst  schnell  ausgedempft  wird.  Haufig  werden 
Fliissigkeitsdampfer  unterschiedlicher  Art  verwendet.  Im  Prinzip  handelt  es  sich  dabei  um 
ein  abgeschlossenes  Gehhuse  in  dem  eine  Fliissigkeit  periodische  Bewegungen  vollfiihren 
kann.  Fiir  EQUATOR-S  wurde  ein  Fliissigkeits-Nutationsdampfer-System  gewahlt  (Bild 
lb)  das  sich  schon  beim  AMPTE-Satelliten  bewahrt  hatte  [3].  Dieser  Nutationsdempfer 
wurde  am  hiesigen  Institut  fiir  Raumfahrttechnik  theoretisch  untersucht  und  gebaut.  Er 
hat  seinen  Anspriichen  voll  geniigt  und  die  nach  dem  EinschuB  des  Satelliten  in  den 
Erdorbit  aufgetretenen  Nutationsschwingungen  innerhalb  weniger  Minuten  ausgedampft. 
Leider  scheiterte  die  Mission  von  EQUATOR-S  spater  aufgrund  von  Latch-up-Effekten  in 
einem  Mikroprozessor  RAM-Baustein. 

Um  eine  gute  Nutationsdampfung  zu  gewahrleisten  muB  bei  der  Auslegung  des 
Fliissigkeits-Nutationsdampfers  Frequenz  und  Dampfung  der  schwingenden  Fliissigkeit  so 
abgestimmt  werden,  dafi  im  gekoppelten  System  die  Nutationsschwingung  moglichst  gut 
gedampft  wird.  Dazu  wird  zunachst  die  Fliissigkeit  im  Dampfer  als  Schwingungssystem 
modelliert.  Gemafi  [4,  5,  7,  8]  geht  man  in  zwei  Schritten  vorg:  Frequenzbestimmung  des 
vereinfachten  reibungsfreien  Modells  und  dann  durch  Einfiihren  eines  Dampfungsterms  die 
Bestimmung  von  gedampfter  Frequenz  und  Dampfung  der  Fliissigkeitsschwingung.  Dieses 
Modell  wird  dann  in  den  nutierenden  Satelliten  integriert  und  die  Nutationsdampfung 
bestimmt.  Man  erhalt  damit  einen  zufriedenstellend  funktionierenden  Nutationsdampfer. 
Die  Funktionstiichtigkeit  eines  solchen  Dampfers  kann  man  auch  im  Erde-Labor  an  einem 
Pendel  testen.  Wir  haben  dazu  ein  luftgelagertes,  also  sehr  reibungsarmes,  Pendel  be- 
nutzt. 

Vergleiche  der  Dampfungswerte  aus  Theorie  und  Auswertungen  der  an  Bord  von 
EQUATOR-S  aufgezeichneten  Nutationsschwingungen  zeigten  jedoch  deutliche  Abwe- 
ichungen.  Dies  hat  uns  veranlaBt,  abweichend  von  [4,  5,  8],  in  einem  dritten  Schritt 
die  tatsachliche  Masse  der  bewegten  viskosen  Fliissigkeit  abzuschatzen;  denn  urspriinglich 
wurde  vereinfachend  vorausgesetzt,  dafi  die  gesamte  reibungsfrei  betrachtete  Fliissigkeit  in 
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dem  betrachteten  Rohr  des  Nutationsdampfers  an  der  schwingenden  Bewegung  teilnimmt. 
Mit  dieser  Betrachtung  erreichen  wir  gut  ubereinstimmende  Ergebnisse  von  Theorie  und 
Messung. 


2  Grundlegende  Gleichungen 


Die  Rotationsbewegung  eines  Satelliten,  der  um  eine  seiner  Haupttragheitsachsen  (z- 
Achse)  spinstabilisiert  ist,  wird  durch  die  Euler- Gleichungen 

Ixx^x  (,^-yy  1 zz)^y^z  Mx  (1) 

lyy^y  i^zz  dxxs)ujz^x  My  (2) 

Izz^z  ( I XX  dyy')LOxU!y  Mz.  (3) 

beschrieben.  Dabei  bezeichnen  Ixx,  Iyy  und  Izz  die  Haupttragheitsmomente  des  Satel¬ 
liten,  ux,  ujy  und  uz  die  Komponenten  der  Winkelgeschwindigkeit  und  Mx,  My  und 
Mz  die  Komponenten  der  aufieren  Momente.  Wirken  keine  auBeren  Momente  auf  den 
Satelliten  und  gilt  bei  Spinstabilisierung  um  die  z-Achse  \ojx\  -C  \uz\,  l^j  <C  \u>z\,  dann 
liefert  die  linearisierte  Gl.  (3),  d.h.  uz  =  0  eine  konstante  2-Komponente  uz  =  uz0  der 
Winkelgeschwindigkeit.  Aus  (1)  und  (2)  folgt  dann  die  Nutationsfrequenz 


Aus  obiger  Nutationsfrequenz  erkennt  man,  dafi  fur  Izz  >  Ixx,  Izz  >  Iyy  Stabilitat  der 
Rotationsbewegung  bereits  im  ungedampften  Fall  gewharleistet  ist.  Die  unerwiinschte 
Nutation  kann  mittels  eingebautem  Dampfungssystem  reduziert  bzw.  zum  Verschwinden 
gebracht  werden. 

Im  folgenden  wird  ein  Nutationsdampfungssystem  mit  zwei  gleichen  symmetrisch  ange- 
brachten  Fliissigkeitsdampfern  betrachtet  (Bild  1).  Jeder  Dampfer  besteht  dabei  aus  zwei 
Rohren  (” “Tubes”)  mit  kreisformigem  Querschnitt  vom  Durchmesser  dT  und  der  Lange 
l  die  durch  zwei  kurze  Rohre  (’’‘Caps”)  mit  ebenfalls  kreisformigem  Querschnitt  vom 
Durchmesser  dc  und  der  Lange  b  miteinander  verbunden  sind.  Die  Dampfer  sind  mit 
Fliissigkeit  der  Dichte  g  und  Viskositat  v  teilweise  gefullt:  jeweils  ein  Tube  vollstandig 
und  die  beiden  Caps  bis  zur  Hohe  b/ 2.  Die  gedampfte  Schwingung  der  Fliissigkeit  im 
Dampfer  wird  beschrieben  durch  [7] 

z  +  2(u>o  z  +  Uq  z  =  0  (5) 


wobei  ro  der  Abstand  des  Dampferschwerpunktes  von  der  Spinachse  des  Satelliten  ist  und 
r0u^0  die  auf  die  Dampferfliissigkeit  wirkende  Beschleunigung  darstellt. 
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Die  Nutationsdampfung  gewinnt  man  dann  aus  der  Euler- Gleichung  fiir  einen  Satelliten 
mit  veranderlichen  Tragheitsmomenten 

(Is  +  lD)  to  +  to  x  (Is  +  I/))  w  +  i^w  +  Kw-l-L  =  Ma.  (7) 

In  (7)  bezeichnen  Is  den  zeitlich  unveranderlichen  Anted  des  Satellitentragheitsmomentes, 
Id  das  zeitlich  veranderliche  Tragheitsmoment,  K  den  Coriolistensor,  L  die  Rela- 


tivbeschleunigung  und  MQ  das  auf  den  Satelliten  wirkende  aufiere  Moment.  Werden 
die  beweglichen  Fliissigkeitsmassen  der  beiden  Dampfer  durch 

ri,2  =  (±f0cosa,  ±r0sina,  h  ±  z)T  (8) 

beschrieben,  wobei  zur  Abkiirzung  f0  =  r0  +  6/2  gesetzt  ist,  dann  berechnet  man  den 
Coriolistensor  aus 

K  =  Ki  +  K2  =  kglAT[{rirJ  -  hif)  +  (r2rj  -  r2rj)]  (9) 

und  fiir  die  Relativbeschleunigung  gilt 

L  =  Li  +  L2  =  kglAT( n  x  r:  +  r2  x  r2),  (10) 


wobei  At  =  kStI^  die  Querschnittsflache  des  Tube  bezeichnet.  Damit  lauten  die  lin- 
earisierten  Bewegungsgleichungen  fiir  das  System  Satellit-Dampfer  [8] 

IiA  —  kglAxfl  sin(2a)  c oy  +  ( Izz  —  Iyy)  tozotoy  +  kglAxfl  sin(2a)  toz0tox  + 

+  2kglATf  o  sin  a(z  +  loIqz)  =  0,  (11) 


lyyUy  -  kglATrl  sin(2a)  tox  +  (Ixx  -  Izz)  toz0tox  -  kglATrl  sin(2a)  toz0toy  - 

-  2 kglATr0  cos  a(z  +  tol0z)  =0  (12) 

und 

z  +  2(u>0z  +  lOqZ  +  [sinacuz  -  cosaay,  +  (cosau;*  +  sin  atoy)  toz0]  r0  =  0.  (13) 

Fiir  die  Aufteilung  der  TYagheitsmomente  in  unveranderliche  (Ixx,  Iyy,  Izz)  und 
veranderliche  (IXXtD,  Iyy,D)  IXZtD)  Anteile  gilt  dabei 


(14) 

(15) 

(16) 

(17) 

(18) 
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ein  Korrekturkaktor  eingefiihrt  wird.  Ohne  diesen  Korrektrufaktor  (18)  wiirde  vorausge- 
setzt,  daft  sich  das  gesamte  Flussigkeitsvolumen  des  Tube  Vt  =  At  l  wie  eine  reibungsfreie 
Fliissigkeit  im  Tube  bewegt.  Wegen  der  Viskositat  der  Flussigkeit  bewegt  sich  jedoch  ein 
geringeres  Flussigkeitsvolumen  Vv  <  VT  im  Tube.  Ein  Vergleich  mit  den  Aufzeichnungen 
an  Bord  von  EQUATOR-S  zeigt,  dab  die  theoretische  Zeitkonstante  mit  Beriicksichtigung 
diese  Korrekturfaktors  gut  mit  der  gemessenen  iibereinstimmt,  ohne  Korrekturfaktor  je¬ 
doch  nur  etwa  halb  so  groB  ist.  Ein  solcher  Korrekturfaktor  wurde  in  fruheren  Arbeiten 
[4,  5,  8]  nicht  beriicksichtigt. 


3  Stabilitatsbetrachtung  und  Nutationsdampfung 

Fiir  die  Stabilitatsuntersuchung  und  zur  Berechnung  der  Nutationsdampfung  liefert 
ein  est-Ansatz  im  obigen  linearen  Differentialgleichungssystem  (11-13)  ein  lineares  alge- 
braisches  Gleichungssystem,  dessen  Determinante  das  charakteristische  Polynom  vierten 
Grades  liefert.  Die  Wurzeln  dieses  Polynoms  liefern  vier  komplexe  Frequenzen  s  —  a±iu, 
deren  Realteil  a  (<  0  bei  Stabilitat)  die  Dampfung  der  Schwingung  angibt.  Aus  dem  Ver- 
schwinden  der  Determinante,  also 

I xx&  T  sin  2 a  ( I zz  lyy)  ^zO  ( Ip  sin  2a)s 

(lxx  -  Izz )  uz0  -  (IF  sin  2a)  s  lyys  -  u)z0IF  sin  2a 

oleosa  +  (sin  a)s  —  (cos  a)s  +  tu2osina 

folgt  das  charakteristische  Polynom 

-f-  G3 4-  0*2 512  -t-  0-iS  -f-  gq  =  0)  (20) 


2 (s2  +  uPz0)IF  sin  a 
—2 (s2  +  uz0 )IF  cos  a 
s2  +  2£u/o  s  +  w  q 


=  0 


(19) 


mit 


a4  =  (1 -/*)(!-/*),  a3  =  2(&(1  —  1*1*), 


&2  =  1  -f-CU2(l  —  I* I  )  + 
a\  =  2CJJJ 

1*7 


3 1*7  -  2(1*  +  7)  +  2(1  +  I)  -  hi  -  I2I 


(h~l)(I2-l) 


1*7 


und 


do  =  u 


1  - 


(h  -  1  )(I3  -  1)J 


+ 


(h  -  1  )(h  ~  1). 

21*7  -  (I*  +  7)  +  2(1  +  7)  -  hi  -  IJ 


(h  -  i)2(/2  -  i)2 


Dabei  wurde  s  =  substituiert  und  zur  Abkiirzung 


<4r  =  (h  ~  1  )(h  -  l)^2o  (~  4),  n;  =  I  =  7~'  I  =  T~ 

lyy  *  xx 

I*  -  2 1  cos2  a,  7  =  27 sin2  a,  h  =  =j7-,  h  =  und  IF  =  gATlrl 

1  xx  yy 

gesetzt.  Die  Nutationsdampfung  wird  durch  den  Realteil  a  der  zur  Nutationsschwin- 
ung  gehorigen  komplexen  Frequenz  s  ausgedriickt.  Eine  anschaulichere  GroBe  fiir  das 
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Abklingverhalten  der  Nutation  ist  die  Zeitkonstante  r  =  -1/a.  Diese  gibt  an  nach 
welcher  Zeit  der  Nutationswinkel  auf  das  1/e  ~  0.37-fache  (rund  1/3)  abgeklungen  ist.  Fiir 
diese  Zeitkonstante  kann  man  mittels  der  ”‘Energy-Sink”-Methode  eine  Naherungsformel 
gewinnen  [4,  8] 


T  = 


_ Ixjjzz  -  Ixx)  [(t»2  ~  wl)2  +  4C 2wjwj] _ 

^C,ojQlzzkQATlflw\[{vji  -  N)2  sin2  a  +  (wxN  -  l)2  cos2  a] 


(21) 


mit 


N  = 


N 


lx  x{lzz  Ixx) 

TJJ7z  -  Iyy)  ’ 


uq  =  — ,  und 
w2o 


re2 


CUo 

^zO 


Es  sei  erwahnt,  da8  fiir  b  <^r0  und  bei  kleiner  Fliissigkeitsmasse  in  den  obigen  Gleichungen 
r0  =  ( ro  +  6/2)  durch  r0  ersetzt  werden  kann. 


4  Test  des  Nutationsdampfers  am  Pendel 


Die  Funktion  des  Nutationsampfer  kann  an  einer  Pendelmaschine  (Bild  31)  getestet  wer¬ 
den.  Es  handelt  sich  dabei  um  ein  physikalisches  Pendel  an  dem  der  Dampfer  in  einem 
wahlbaren  Abstand  von  der  Drehachse  des  Pendels  montiert  werden  kann.  AuCerdem 
kdnnen  Frequenz  und  Energie  des  Pendels  durch  verschiebbare  starre  Massen  eingestellt 
werden.  Fiir  dieses  Pendel  leiten  wir  im  folgenden  die  Bewegungsgleichungen  mit  Hilfe 
der  Lagrangeschen  Gleichungen 


d_  (dT_\  dV_  dD 
dt  \dqn)  +  dqn  +  dqn 


(22) 


her.  Mit  den  zwei  generalisierten  Koordinaten  cp  (Drehwinkel  des  Pendels)  und  z  (Be- 
wegungsstrecke  der  Fliissigkeit  im  Tube)  hat  das  System  Pendel-Dampfer  die  kinetische 
Energie 

T  =  ~Is^  +  ~mD(z  +  Topf  +  ^/o02,  (23) 

die  potentielle  Energie 

V  =  ~rn.£)CUgZ2  +  mDgzshup  +  mDgr0(  1  -  cosy?)  +  Mgrs{  1  -  cosy?)  (24) 
und  die  zugehorige  Dissipationsfunktion  lautet 


D  =  ^mD2(u0z2.  (25) 

Dabei  bezeichnen  Is  das  Tragheitsmoment  der  starren  Massen,  mo  =  kglAr  die 
Fliissigkeitsmasse  im  Tube,  r0  den  Abstand  der  Fliissigkeit  von  der  Drehachse,  I0  = 
Y^mDl2  das  Tragheitsmoment  der  als  starr  angenommenen  Fliissigkeitssaule  um  ihre 
zur  Saule  senkrechten  Schwerpunktsachse,  cj0  die  ungekoppelte  Eigenfrequenz  der 
Fliissigkeitssaule,  g  die  Schwerebeschleunigung,  M  den  starren  Anteil  der  Pendelmasse 
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und  rs  den  Schwerpunktsabstand  der  starren  Pendelmasse  vom  Drehpunkt.  Hier  wurde 
wieder  fur  die  bewegte  Masse  mo  ein  entsprechender  Korrekturfaktor 

k  =  16W^I  (26) 

gdT 

eingefiihrt. 

Anwendung  des  Lagrangeschen  Formalismus  und  Linearisierung  liefert  die  beiden  Bewe- 
gungsgleichungen  (Q1j2  =  0) 

Ig<P  +  MGrsgv  +  mDfQZ  +  mDgz  =  0  (27) 

und 

z  +  2(uqz  +  UqZ  +  r0<p  +  gp  =  0,  (28) 


1  +  7 


und 


wobei  IG  =  Is  +  Id  das  Tragheitsmoment  des  Gesamtsystems  um  die  Drehachese  des 
Pendels,  Id  =  Io+morl  das  Tragheitsmoment  der  Fliissigkeit  um  die  gleiche  Achse,  MG  — 
M  +  mo  die  Gesamtmasse  des  Pendels,  rs  den  Schwerpunktsabstand  des  Gesamtsystems 
von  der  Drehachse  des  Pendels  und  g  die  Schwerebeschleunigung  bezeichnen.  Fiir  letzteren 
gilt  die  Beziehung  MGfs  =  Mrs  +  ra^Fo-  Mit  der  Pendelfrequenz 

/  MGrsg 

"» =  y  “jy  <30) 

und  der  Substitution  s  =  upS  erhalten  wir  mit  dem  est  =  e^^-Anstatz  das  charakteris- 
tische  Polynom 


mDrl 


54  +  2C^53+(l  +  l-^^]52  +  2C^+^-^  =  0.  (31) 


ul  2 mDr0\  o2 


ul  mDIG 


V  Ig  J~  V*  '  ul  MGrs  J  '  '  uj  Mir \  . ' 

Die  Zeitkonstante  r  erhalt  man  aus  der  zum  Pendel  gehorigen  Losung  S  des  obigen  eharak- 
teristischen  Polynoms  (31)  als 

1  1 


u;pRe  S 


5  Numerische  Auswertung 

5.1  Dampfung  der  Nutationsschwingung  des  Satelliten 

Wie  schnell  eine  Nutation  durch  den  Dampfer  abgebaut  wird  kann  durch  die  oben 
angegebene  Zeitkonstante  r  ausgedruckt  werden.  Diese  Zeitkonstante  gibt  an  nach  welcher 
Zeit  der  Nutationswinkel  auf  das  1/e-fache,  also  etwa  den  dritten  Teil  abgeklungen  ist. 
Um  r  zu  bestimmen,  mussen  die  Wurzeln  des  charakteristischen  Polynoms  (20)  numerisch 
berechnet  werden.  Aus  dem  zur  Nutationsschwingung  gehorigen  Losungspaar  wurde  die 
Zeitkonstante  r  =  — 1/Re(s)  ermittelt.  Um  den  EinfluB  der  verschiedenen  Parameter  zu 
erkennen,  wurde  die  Zeitkonstante  als  Funktion  eines  der  Parameter  aufgetragen. 
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5.1.1  Maximale  Dampfung 

Zur  Ermittlung  der  starksten  Dampfung  kann  zunachst  einmal  -  mit  eventuellen  Ein- 
schrankungen  fiir  die  Parameter  -  der  Minimalwert  der  Zeitkonstanten  r  berechnet  wer- 
den.  Ein  solches  Minimum  wurde  hier  ermittelt;  allerdings  fur  den  numerisch  einfacher 
zu  bestimmenden  Naherungswert  f .  Hierzu  haben  wir  folgendes  angenommen:  D”ampfer 
aus  Aluminium  (Dichte  2700  kg/m3)  mit  Rohrwandstarken  von  1.5  mm  fiir  den  Tube  und 
2  mm  fiir  das  Cap,  einen  Satelliten  mit  den  Tragheitsmomenten 

Ixx  =  33.34  kgm2,  Iyy  =  34.80  kg  m2,  Izz  =  45.35  kgm2, 

Spinfrequenz,  Einbauhohe  sowie  Einbauwinkel  der  Dampfer 

cu2o  =  48  rpm,  H  —  0,  a  =  0, 

Dichte  und  kinematische  Zahigkeit  der  Dampferfliissigkeit 

q  =  1841  kg/m3,  v  =  1.3  cSt. 

Mit  den  Einschrankungen  fiir  die  restlichen  geometrischen  Parameter 

40  mm  <  dc  <  60  mm  fiir  den  Innendurchmesser  der  Cap-Rohre, 

40  mm  <  b  <  80  mm  fiir  die  Lange  der  Cap-Rohre, 

8  mm  <  dr  <  10  mm  fiir  den  Innendurchmesser  der  Tube-Rohre, 

400  mm  <  l  <  500  mm  fiir  die  Lange  der  Tube-Rohre  und 

400  mm  <  r0  <  565  mm  fiir  den  Einbauradius  der  Dampfer. 

liefert  die  Optimierung  das  Ergebnis 

Durchmesser  Cap  dc  =  46.08  mm 

Lange  des  Cap  b  =  80.00  mm  (Einschrankungsgrenze) 

Durchmesser  Tube  dx  =  10.00  mm  (Einschrankungsgrenze) 

Lange  des  Tube  l  —  500.00  mm  (Einschrankungsgrenze) 

Einbauradius  ro  =  565.00  mm  (Einschrankungsgrenze) 

und  als  Minimum  der  Zeitkonstanten 

fmin  =  50.33  Sekunden. 

5.1.2  Einflufi  der  Parameter  auf  die  Nutationsdampfung 

Der  Nutationsdampfer  fiir  den  Satellit  EQUATROR-S  sollte  fiir  Spinraten  bis  50  rpm  und 
verschiedene  Tragheitsmomente  (Ausklappen  der  Booms  nach  dem  Aufspinnen)  wirken. 
Der  Dampfer  wurde  daher  wie  folgt  ausgelegt 

Durchmesser  Cap  dc  =  45  mm 

Lange  des  Cap  b  =  80  mm 

Durchmesser  Tube  dT  =  9  mm 

Lange  des  Tube  l  =  500  mm 
Einbauradius  r0  =  565  mm. 
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Die  restlichen  Parameter  sind  jene  bereits  oben  angegebenen.  In  den  Bildern  2  und  3 
wurde  der  Korrekturfaktor  k  in  Abhangigkeit  des  Tube-Durchmessers  dr  (Bild  2)  bzw. 
des  Cap-Durchmessers  dc  (Bild  3)  aufgetragen.  Aus  Bild  3  kann  man  entnehmen,  dab 
der  Korrekturfaktor  fur  die  oben  gewahlten  Parameter  und  dr  —  9  mm  den  Wert  k  ~  1/3 
hat,  d.h.  dab  die  bewegte  Fliissigkeitsmasse  nur  ein  Drittel  der  im  Tube  befindlichen 
Fliissigkeitsmasse  qAtI  betragt.  Weiter  kann  man  aus  Bild  2  erkennen  dab  der  Korrektur¬ 
faktor  mit  wachsendem  Tubedurchmesser  weiter  abnimmt  und  fur  kleinere  d?  grober  wird 
(k  =  0.6  fiir  dr  =  8  mm.  Die  Abhangigkeit  von  k  vom  Cap-Durchmesser  ist  umgekehrt, 
d.h.  mit  wachsendem  dc  wachst  auch  der  Korrekturfaktor  (k  ~  0.65  fiir  dc  =  60 mm). 
Wie  sich  die  Korrektur  auf  die  Zeitkonstante  auswirkt  ist  in  Bild  4  gezeigt.  Hier  sind  fiir 

variablen  Tubedurchmesser  8  mm  <  dT  <  9.5 mm  die  korrigierte  Zeitkonstante  r  ( - ), 

die  korrigierte  Naherungslosung  f  ( - ),  die  unkorrigierte  (ohne  Korrekturfaktor  k,  d.h. 

k  =  1)  Zeitkonstante  r(ohne  k)  ( - )  und  die  unkorrigiert  Naherungslosung  f  (ohne  k) 

(••••)  sowie  das  Ergebnis  aus  der  Anbordaufzeichnung  von  EQUATOR-S  texv  aufgetragen. 
Man  erkennt  den  deutlichen  Unterschied  zwischen  den  korrigierten  und  unkorrigierten 
Werten  fiir  die  Zeitkonstante.  Fiir  die  gewahlten  Parameter  dy  =  9  mm  ist  die  korrigierte 
Zeitkonstante  etwa  doppelt  so  grob  wie  die  unkorrigierte.  Der  gemessene  Wert  stimmt 
dabei  sehr  gut  mit  der  korrigierten  Zeitkonstanten  r  iiberein.  Dies  zeigt,  dab  fiir  eine  quan- 
titativ  gute  Voraussage  der  Korrekturfaktor  beriicksichtigt  werden  mub.  In  den  Bildern 

5-16  ist  die  Zeitkonstante  r  ( — )  und  deren  Naherungswert  t  ( - )  als  Funktion  jeweils 

eines  der  Parameter  aufgetragen.  Der  gemessene  Wert  an  Bord  von  EQUATOR-S  ist  eben- 
falls  wieder  mit  eingezeichnet.  In  Bild  5  wurde  das  Tragheitsmoment  um  die  x-Achse  Ixx 
des  Satelliten  variert.  Hier  kann  man  erkennen,  dab  Anderungen  des  Tragkeitsmomentes 
die  Zeitkonstante  deutlich  beeinflussen.  Verringert  sich  Ixx  so  steigt  die  Zeitkonstante 
an;  eine  Anderung  von  weniger  als  10%  des  Tragheitsmomentes  kann  bereits  eine  Ver- 
doppelung  der  Zeitkonstante  bedeuten.  Ein  anliches  Verhalten  gilt  auch  bei  Anderung 
des  Tragheitsmoment  um  die  y-Achse  Iyy  (Bild  6).  Es  zeigt  sich  bereits  hier,  dab  die 
endgultigen  Tragheitsmomente  des  Satelliten  von  vornherein  bekannt  sein  miibten,  wollte 
man  die  bestmogliche  Nutationsdampfung  erzielen.  Der  Einflub  des  Tragheitsmomentes 
Izz  um  die  Spinachse  ist  hier  umgekehrt,  d.h.  mit  kleiner  werdendem  Izz  nimmt  auch  die 
Zeitkonstante  r  ab  (Bild  7)  und  hat  fiir  Izz  «  44 kgm?  einen  Minimalwert  von  r  «  45s 
(Bilder  7  und  8).  In  diesen,  wie  auch  in  den  folgenden  Bildern  erkennt  man  auch,  dab  der 
Naherungswert  t  fiir  die  Zeitkonstante  r  recht  gut  mit  dieser  iibereinstimmt  und  daher 
fiir  die  Auslegung  eines  Nutationsdampfers  geeignet  ist. 

Der  Einflub  der  Spinrate  ist  in  Bild  9  dargestellt.  Mit  kleiner  werdender  Spinrate  nimmt 
die  Zeitkonstante  zu  und  hat  fiir  uz o  =  30  rpm  den  Wert  r  =  95  s.  Die  Nutationsdampfung 
hangt  auch  stark  vom  Abstand  r0  des  Dampfers  von  der  Spinachse  ab.  Dies  kann  aus 
Bild  10  entnommen  werden.  Wie  zu  erwarten  steigt  die  Zeitkonstante  mit  abnehmendem 
r0  und  hat  fiir  r0  =  0.45  m  den  Wert  r  =  150  s,  ist  also  mehr  als  doppelt  so  grob  wie 
fiir  den  gewahlten  Abstand  r0  =  0.565  m  mit  einem  r  =  72  s.  Einen  starken  Einflub 
auf  die  Dampfung  hat  auch  die  Lange  des  Tube.  In  Bild  11  ist  die  Zeitkonstante  r  als 
Funktion  der  Tubelange  l  aufgetragen.  Bei  dem  hier  gewahlten  festen  Parametersatz  zeigt 
sich  interessanterweise,  dab  mit  kleiner  werdender  Tubelange  die  Zeitkonstante  fallt,  eine 
kiirzere  Tubelange  also  giinstiger  ware.  Bei  ausgeklapptem  Booms  jedoch,  also  anderen 
Tragheitsmomenten,  verhalt  sich  dies  jedoch  umgekehrt  (siehe  Bild  24).  Der  Einflub  des 
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Tubedurchmessers  dr  wird  in  Bild  12  gezeigt.  Man  sieht,  daB  mit  wachsendem  d?  die 
Zeitkonstante  abnimmt.  Relativ  gering  ist  der  EinfluB  der  Caplange  wie  man  in  Bild 
13  sehen  kann.  Mit  wachsender  Caplange  b  nimmt  die  Zeitkonstante  geringfiigig  ab. 
Der  EinfluB  des  Capdurchmessers  ist  hingegen  starker.  Bild  14  zeigt,  daB  mit  geringer 
werdendem  Capdurchmesser  dc  die  Zeitkonstante  zunachst  fallt  und  fur  dc  ~  41.5  cm  den 
Minimalwert  r  «  60s  erreicht  und  dann  wieder  ansteigt.  Ebenso  steigt  die  Zeitkonstante 
fiir  dc  >  45  cm  und  hat  beispielsweise  fiir  dc  =  50  cm  bereits  den  Wert  r  =  100  s. 

Einenen  entscheidenden  EinfluB  auf  die  Nutationsdampfung  hat  die  Wahl  der  Fliissigkeit. 
Von  Dichte  und  Viskositat  der  Fliissigkeit  hangt  die  Zeitkonstante  sehr  stark  ab.  Beson- 
ders  drastisch  wirkt  sich  die  Viskositat  auf  die  Zeitkonstante  aus.  Es  sollte  daher  die 
Temperaturabhangigkeit  der  Viskositat  mit  beriicksichtigt  werden.  Bei  EQUATOR-S 
betragt  die  Temperatur  fiir  die  Dampferfliissigkeit  12.7°  C.  Als  Fliissigkeit  wurde  fiir  den 
EQUATOR-S-Nutationsdampfer  PP3  gewahlt.  Bild  15  zeigt  wie  die  Zeitkonstante  von 
der  Fliissigkeitsdichte  g  abhangt.  Fiir  eine  Fliissigkeit,  deren  Dichte  nur  etwa  halb  so  groB 
wie  die  gewhalte  (g  =  1841  kg/m3)  ware,  also  etwa  die  Dichte  von  Wasser  hatte,  ware  die 
Zeitkonstante  etwa  doppelt  so  groB.  Bild  16  zeigt,  daB  der  EinfluB  der  kinematischen 
Viskositat  v  noch  starker  ist.  Eine  geringfiigig  kleinere  Viskositat  v,  etwa  durch  hohere 
Temperatur  der  Fliissigkeit,  bewirkt  bereits  eine  deutlich  hohere  Zeitkonstante  r. 

Die  Bilder  17-29  zeigen  die  Ergebnisse  fiir  EQUATOR-S  mit  ausgeklappten  Booms.  Die 
Tragheitsmomente  des  Satelliten  sind  dann 

Ixx  =  38.00  kgm2,  Iyy  =  39.37  kgm 2 ,  Izz  =  54.67  kgm2. 

Die  restlichen  Parameter  sind  jene  von  oben.  In  den  Bildern  17  und  18  ist  wieder  der 
Korrekturfaktor  k  in  Abhangigkeit  des  Tubedurchmessers  dx  bzw.  Capdurchmessers  dc 
gegeben.  Das  Verhalten  ist  wieder  ahnlich  jenem  in  den  Bildern  2  und  3.  Der  Kor¬ 
rekturfaktor  hat  jetzt  aufgrund  der  etwas  hoheren  Nutationsfrequenz  einen  geringfiigig 
hoheren  Wert.  Die  Zeitkonstante  ist  jetzt  mit  r  =  345  s  bei  ausgeklappten  Booms 
wesentlich  groBer.  Das  Verhalten  bei  geringer  Anderungen  der  Haupttragheitsmomente 
ist  jedoch  wieder  ahnlich  wie  oben  (vgl.  Bilder  19  bis  21  mit  den  Bildern  5  bis  7).  Die 
Anderung  der  Zeitkonstanten  bei  geringer  Anderung  der  Spinfrequenz  ujzq  verhalt  sich 
jetzt  umgekehrt,  d.h.  mit  kleiner  werdendem  u>zq  wird  auch  die  Zeitkonstante  kleiner 
(Bild  22).  Bei  r  »  25  rpm  wurde  r  minimal  und  dann  wieder  ansteigen  (hier  nicht 
gezeigt).  Wie  mit  eingeklappten  Booms  verhalt  sich  die  Abhangigkeit  der  Zeitkonstan¬ 
ten  vom  Dampferabstand  r0  auch  hier  wieder  (Bild  23).  In  Bild  24  kann  man  sehen, 
daB  fiir  die  gewahlte  Tubelange  l  =  0.5  m  die  Zeitkonstante  fiir  EQUATOR-S  mit  aus¬ 
geklappten  Booms  nahezu  ihr  Maximum  erreicht,  d.h.  ein  kiirzerer  oder  langerer  Tube 
wahre  jetzt  giinstiger.  Gleiches  gilt  auch  fiir  den  Tubedurchmesser  (Bild  25)  wahrend 
sich  die  Caplangenanderung  wieder  wie  friiher  auf  die  Zeitkonstante  auswirkt  (Bild  26). 
Anderungen  des  Capdurchmessers  wirken  sich  ebenfalls  ahnlich  wie  oben  aus;  nur  wird 
jetzt  t  nicht  schon  fiir  dc  im  dargestellten  Bereich  40  mm  <  dc  <  50  mm  minimal  (Bild 
27).  Die  Abhangigkeit  der  Nutationsdampfung  von  den  Fliissigkeitsparametern  Dichte  g 
und  Viskositat  v  ist  wieder  wie  friiher  (vgl.  Bilder  28  und  29  mit  den  Bildern  15  und  16). 

Wie  die  obigen  Ergebnisse  zeigen,  war  der  Nutationsdampfer  nicht  so  ausgelegt  daB  die 
bestmogliche  Dampfung,  d.h.  kleinste  Zeitkonstante  r  erzielt  wurde.  Das  liegt  daran,  dafi 
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die  Tragheitsmomente  von  EQUATOR-S  nicht  von  Anfang  an  fest  vorgegeben  werden  kon- 
nten,  daB  auBerdem  eine  gute  Dampfung  sowohl  mit  eingeklappten  als  auch  ausgeklappten 
Booms  gegeben  sein  mufite  und  daB  schlieBlich  gewisse  Einschrankungen  gegeben  waren, 
etwa  bei  der  Wahl  von  Tube-  und  Capdurchmesser.  Wie  aus  den  Aufzeichnungen  an  Bord 
von  EQUATOR-S  zu  sehen  ist  (siehe  unten),  hat  der  Nutationsdampfer  jedoch  seine  Auf- 
gabe  gut  erfiillt. 

5.1.3  Aufzeichnung  der  Nutationsdamfung  an  Bord  von  EQUATOR-S 

Bild  30  zeigt  das  Abklingen  der  Nutationsamplitude  von  EQUATOR-S  zu  Beginn  der 
Mission.  Der  Anfangswinkel  der  Nutation  betrug  3.84°.  Aus  der  Aufzeichnung  erhalt 
man  fiir  die  Zeitkonstante  den  (gemittelten)  Wert  texp  =  (72  ±  6)s.  Aus  der  Rechnung 
folgt  fiir  die  Parameter 

Ixx  —  33Mkgm2,  Iyy  =  34.80  kgm2,  Izz  =  45.35  kgm2 

lozo  =  48  rpm,  ro  —  0.565  m  H  =  0,  a  =  0 
dT  =  9mm,  l  =  39  cm,  dc  =  43  mm,  b  =  8  cm 
g  =  1841  kg/m3,  u=  1.3  cSt 

eine  Zeitkonstante  r  =  72.3  s.Das  ist  eine  sehr  gute  Ubereinstimmung  mit  dem  gemessenen 
Wert. 

5.2  Dampfung  der  Pendelschwingung  mittels  Nutationsdampfer 

Der  Nutationsdampfer  kann  im  Erdlabor  an  einem  Pendel  getestet  werden.  Die  Zeitkon¬ 
stante  laBt  sich  hier  sowohl  numerisch  aus  der  charakteristischen  Gleichung  (31)  wie 
auch  experimentell  ermitteln.  Ein  Vergleich  beider  Ergebnisse  erlaubt  Aussagen  iiber 
die  Gute  des  mathematischen  Modells.  Das  verwendete  Pendel  ist  luftgelagert  und 
hat  somit  geringe  Reibung.  Das  Abklingen  der  Pendelschwingung  kommt  daher  nahezu 
ausschlieBlich  durch  den  Dampfer  zustande.  Am  Pendel  konnen  zwei  starre  Massen  an 
zwei  Aluminiumrohren  oberhalb  und  unterhalb  des  Pendeldrehpunktes  verschiebbar  ange- 
bracht  werden  (Bild  31).  Dadurch  lassen  sich  Schwerpunktsabstand  und  Tragheitsmoment 
und  damit  Pendelfrequenz  und  -energie  variieren.  Die  Daten  fiir  dieses  Pendel  konnen  aus 
Bild  32  entnommen  werden. 

5.2.1  Numerische  Ergebnisse  der  Dampfung 

Fiir  die  Versuche  am  Pendel  wurde  ein  Dampfer  aus  Plexiglas  mit  den  Abmessungen 
dT  =  8  mm,  l  —  46  cm,  dc  —  34  mm,  und  b  —  11.25  cm 
verwendet.  Als  Fliissigkeit  wurde  Wasser  benutzt  und  hierfiir 

g  =  1000  kg/m3,  und  u  =  1  cSt 
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angenommen.  Das  Dampergehause  hatte  eine  Masse  mDG  =  0.4829%  und  ein 
Tragheitsmoment  um  seine  Schwerpunktsachse  (parallel  zur  Drehachse  des  Pendels) 
Idg  =  1.592-10-2  kg  m2.  Mit  den  Abstanden  der  starren  Massen  Lx  =  65  cm  (oberhalb  des 
Drehpunktes)  und  L2  =  90  cm  (unterhalb  des  Drehpunktes)  erhalt  man  fur  einen  Dreh- 
punktsabstand  r0  =  67  cm  des  Dampferschwerpunktes  eine  theoretische  Zeitkonstante 
(aus  Gl.  (31)  ) 

rth  =  38.7  s. 

Das  Experiment  lieferte  fur  diesen  Fall  (siehe  unten) 

rei  =  (38  ±  3)  s. 

Damit  hat  man  eine  gute  Ubereinstimmung  von  Theorie  und  Experiment.  Audi  hier 
zeigte  sich  wieder  die  Notwendigkeit  des  Korrekturfaktors.  Ohne  diesen  wiirde  die  Rech- 
nung  eine  Zeitkonstante  von  nur  11.2  s  liefern,  also  erheblich  vom  experimentellen  Wert 
abweichen.  Durch  den  Korrekturfaktor  k  wurde  eben  beriicksichtigt,  dafi  fiir  die  be- 
wegte  Fliissigkeitsmasse  nicht  das  gesamte  Fliisssigkeitsvolumen  des  Tube  gesetzt  werden 
darf,  sondern  aufgrund  der  viskosen  Rohrstromung  ein  zu  berechnender  geringerer  Wert 
genommen  werden  mufi.  Dabei  mufi  jedoch  beachtet  werden,  dafi  nach  den  Voraussetzun- 
gen  der  Theorie  die  Parameter  nicht  vollig  beliebig  gewahlt  werden  konnen,  sondern  daB 
beispielsweise  das  Verhaltnis  dr/dc  klein  und  auch  laminare  Stromung  gegeben  sein  mufi. 

5.2.2  Experimentelle  Ergebnisse 

Fiir  das  Experiment  wurde  der  oben  beschriebene  Dampfer  verwendet.  In  den  Bildern 
33  und  34  sind  die  Pendelschwingungen  als  Funktion  der  Zeit  fiir  das  Pendel  mit 
Dampfergehause  ohne  Fliissigkeit  (Bild  33)  und  mit  wassergefiilltem  Dampfer  (Bild  34) 
aufgezeichnet.  In  Bild  33  ist  die  Dampfung  aufgrund  der  Luftlagerung  des  Pendels  auBerst 
gering.  Die  starke  Dampfung  der  Schwingung  in  Bild  34  wird  also  fast  ausschliefilich  durch 
den  Fliissigkeitsdampfer  erzielt.  Eine  Auswertung  der  Aufzeichnung  bringt  fiir  kleine 
Schwingungsamplituden  eine  Zeitkonstante  texp  =  (38  ±  3)  s  und  stimmt  gut  mit  dem 
Wert  rth  =  38.7s  der  Theorie  iiberein  (siehe  oben).  Fiir  den  Vergleich  mit  der  linearen 
Theorie  ist  dabei  darauf  zu  achten,  daB  im  Experiment  kleine  Schwingungsamplituden 
gegeben  sind. 


Anhang 

A  Herleitung  des  Korrekturfaktors 

Zur  Herleitung  des  Korrekturfaktors  k  in  (18)  und  (26)  vergleichen  wir  das  Durchflufivol- 
umen  der  reibungsfreien  Stromung  im  Rohr  mit  jenem  der  viskosen  Rohrstromung.  Es 
soil  also  ermittelt  werden,  wie  mittels  eines  Korrekturfaktors  k  das  Durchstromvolumen 
einer  reibungsfreien  Fliissigkeit  korrigiert  werden  mufi,  so  daB  gilt 

Vtv  =  kVrrf ,  (33) 
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wenn  VTrf  dem  Volumen  der  reibungsfreien  Fliissigkeit  entsprache  und  VTv  das  Volu- 
men  der  viskosen  Fliissigkeit  bezeichnet.  Dem  DurchfluB  einer  reibungsfreien  Fliissigkeit 
entsprache  eine  Amplitude  (rf  der  Fliissigkeitsoberflache  im  Cap  und  jenem  einer  viskosen 
Fliissigkeit  eine  Amplitude  (v.  Es  gilt  also  wegen  (vAc  =  k(rjAc 

C»  =  Kr /•  (34) 

Das  DurchfluBvolumen  der  reibungsfreien  Fliissigkeit  wahrend  einer  viertel  Periode  ware 
nun 

Vrf  =  C rfAc  (35) 

mit  einem  unbekannten  £r/-  Um  dieses  zu  eliminieren  berechnen  wir  fiir  das  gleiche 
Volumen  Vv  =  Vrf  den  DurchfluB  im  viskosen  Fall.  Dazu  nehmen  wir  vereinfacht  zu  jedem 
Zeitpunkt  t  eine  stationare  Rohrstromung  an  und  integrieren  dann  fiber  das  Zeitintervall 
einer  viertel  Periodenschwingung. 

Fiir  die  stationare  Rohrstromung  gilt  die  Geschwindigkeitsverteilung 

w  =  -  R2),  (36) 

ir]l 

wobei  po  und  pi  den  Druck  an  den  Rohrenden,  p  die  dynamische  Vikositat,  l  die  Rohrlange 
(Tubelange) ,  r  die  radiale  Rohrkoordinate  und  R  =  dT/ 2  den  Rohrradius  bezeichnen.  Mit 
den  Auslenkungen  ±(v  cos  (fit)  iiber  der  Gleichgewichtslage  gilt  fiir  den  um  eine  viertel 
Periode  verschobenen  treibenden  Druck  der  Rohrstromung 

Po  =  eg[b/  2  +  C«  sin  (fit)],  pi  =  £>^[6/2  -  sin(fit)]. 

Dabei  steht  Q  fiir  die  Nutationsfrequenz  ojn  beim  Nutationsdampfer  im  Satelliten  und 
fiir  die  Pendelfrequenz  up  beim  Nutationsdampfer  am  Pendel.  Entsprechend  ist  fiir  g  die 
Erdschwere  beim  Pendel  und  tqujzq  beim  Satelliten  zu  setzen.  Fiir  die  Geschwindigkeit 
der  viskosen  Rohrstromung  gilt  dann 


w=M^(r2-JR2)sin(m).  (37) 

Apt 

Das  Durchflufivolumen  der  viskosen  Fliissigkeit  fiir  eine  viertel  Periodenschwingung  ist 
damit  gegeben  durch 


/•7r/2n  fR 

14  =  27 r  /  /  wrdrdt 

J  t=0  J r= 0 


ng(vgR 4 
4plQ 


(38) 


oder  wegen  (34) 


Direkter  Vergleich  mit  (34)  liefert 


v  7TQk(rfgRA 

Apin 


C rfAc  —  Vrf  Vv 


TTQkCrfgR 4 

4pin 


(39) 
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und  damit  wegen  v  —  rjfg,  R  =  dr/2,  Ac  =  7r d2c /4  den  Korrekturfaktor 


1  QulVtd2c 

gdj 


(40) 


mit  Q  =  un,  g  =  tqlo2z0  flir  den  Dampfer  im  Satelliten  und  Q  =  ujp,  g  =  Erdbeschleunigung 
fiir  den  Dampfer  am  Pendel. 
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Nomenklatur 


At  Querschnittsflache  des  Tube 

b  Caplange 

dc  Capdurchmesser 

dr  Tubedurchmesser 

g  Erdbeschleunigung 

h  Abstand  des  Dampferschwerpunktes  von  der  x  -  y-Ebene 

/0  Tragheitsmoment  der  Fliissigkeitssaule  im  Tube  beim  Pendel 

ID  Tragheitsmoment  der  Fliissigkeit  am  Pendel  (Id  =  h  + 

IG  Tragheitsmoment  des  Pendels  (Ig  —  Is  +  Id) 

Is  Tragheitsmoment  der  starren  Masse  des  Pendels 

I5  zeitlich  unveranderlicher  Anted  des  Satellitentragheitsmomentes 

ID  zeitlich  veranderlicher  Anteil  des  Satellitentragheitsmomentes 

I xx,  lyy )  I zz  Haupttragheitsmomente  des  Satelliten  (ohne  schwingende  Dampferfliissigkeit) 
~Ixx,  lyy ,  Izz  Haupttragheitsmomente  des  Satelliten  mit  schwingender  Dampferfliissigkeit 
Ixx  D,  Iyy,D,  zeitlich  veranderliche  Anteile  der  Haupttragheitsmomente 
Izz,D  des  Satelliten  (schwingende  Dampferfliissigkeit) 

k  Korrekturfaktor 

K  Coriolistensor 

L  Relativbeschleunigung 

l  Tubelange 

M  starre  Pendelmasse 

Mq  Gesamtmasse  des  Pendels  (Mg  =  M  +  mo) 

Ma  aufleres  Moment 

Mx,  My,  Mz  Komponenten  der  aufleren  Momente 

r0  Abstand  des  Dampferschwerpunktes  von  der  Spinachse 

r0  Abstand  der  Tubefliissigkeit  von  der  Spinachse  (Fo  =  ro  +  b/2) 

rs  Schwerpunktsabstand  beim  Pendel 


a 

c 

v 

Q 

UN 

Uo 

(D  x  >  ^  z 

<^z0 

UJp 


Winkel  der  Dampferebene  mit  der  z-Achse 

Dampfungsfaktor  der  ungekoppelten  Dampferschwingung 

kinematische  Viskositat  der  Dampferfliissigkeit 

Dichte  der  Dampferfliissigkeit 

Nutationsfrequenz 

ungekoppelte  Dampferfrequenz 

Winkelgeschwindigkeitskomponenten 

Spinfrequenz 

Pendelfrequenz 
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18 


Bild  la  B.  Hausler/  W.  Eidel 


Bild  lb  B.  Hausler/  W.  Eidel 
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Bild  lc  B.  Hausler/  W.  Eidel 


1^  =  33.34  kg  m‘ 
I  =  34.80  kg  m5 
l=  =  45.35  kg  m‘ 
to20  =  5.02655  s'1 

v  =1.3  cSt 
I  =  0.5  m 
rQ  =  0.565  m 

dc  =  45  mm 


•  • 


Zeitkonstante 


Bild  6  B.  Hausler/  W.  Eidel 


Zeitkonstante 


lzz  [kg  m2] 


Bild  8  B.  Hausler/  W.  Eidel 


Zeitkonstante 


I  =  33.34  kg  m2 

XX 

I  =  34.80  kg  m2 
yy  o 

I  =  45.35  kg  m 

r0  =  0.565  m 

I  =  0.5  m 
dT  =  9  mm 

b  =  8  cm 
cL  =45  mm 


_ I _ ' - 1— 

45  50 

coz0  [rpm] 
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Bild  9  B.  Hausler/  W.  Eidel 


Bild  10  B.  Hausler/  W.  Eidel 


Bild  1 1  B.  Hausler/  W.  Eidel 
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Bild  13  B.  Hausler/  W.  Eidel 


Zeitkonstante 


Bild  14  B.  Hausler/  W.  Eidel 


Korrekturfaktor  k 


Bild  18  B.  Hausler/  W.  Eidel 


Zeitkonstante 


I  =  39.37  kg  m2 

yy  2 

la  =54.67  kg  m2 

“zo  =  48  rpm 
rQ  =  0.565  m 

I  =  0.5  m 
dT  =  9  mm 

b  =  8  cm 


i 
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Bild  20  B.  Hausler/  W.  Eidel 


Zeitkonstante 


Bild  21  B.  Hausler/  W.  Eidel 


Bild  22  B.  Hausler/  W.  Eidel 


Zeitkonstante 


dy  [mm] 


Bild  25  B.  Hausler/  W.  Eidel 


Bild  26  B.  Hausler/  W.  Eidel 


800 


Bild  29  B.  Hausler/  W.  Eidel 


Nutationswinkel 

(Aufzeichnung  an  Bord  von  EQUATOR-S) 
_  .  S0D6D2  (Solid); 


Bild  30  B.  Hausler/  W.  Eidel 


Bild  31a  B.  Hausler/  W.  Eidel 


